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ABSTRACT 

We selected a sample of 33 Gamma-Ray Bursts (GRBs) detected by Swift, with known red- 
shift and optical extinction at the host frame. For these, we constructed the de-absorbed and 
K-corrected X-ray and optical rest frame light curves. These are modelled as the sum of two 
components: emission from the forward shock due to the interaction of a fireball with the 
circum-burst medium and an additional component, treated in a completely phenomenolog- 
ical way. The latter can be identified, among other possibilities, as "late prompt" emission 
produced by a long lived central engine with mechanisms similar to those responsible for the 
production of the "standard" early prompt radiation. Apart from flares or re-brightenings, that 
we do not model, we find a good agreement with the data, despite of their complexity and di- 
versity. Although based in part on a phenomenological model with a relatively large number 
of free parameters, we believe that our findings are a first step towards the construction of 
a more physical scenario. Our approach allows us to interpret the behaviour of the optical 
and X-ray afterglows in a coherent way, by a relatively simple scenario. Within this context 
it is possible to explain why sometimes no jet break is observed; why, even if a jet break is 
observed, it is often chromatic; why the steepening after the jet break time is often shallower 
than predicted. Finally, the decay slope of the late prompt emission after the shallow phase is 
found to be remarkably similar to the time profile expected by the accretion rate of fall-back 
material (i.e. oc t~ 5 ^ 3 ), suggesting that this can be the reason why the central engine can be 
active for a long time. 
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The GRB X-ray light curves, as observed by Swift, have shown a 
complexity unforeseen before. Besides the behaviour as observed 
by BeppoSAX after several hours from the trigger, a significant 
fraction of GRBs shows a steep flux decay soon after the end of 
the prompt as seen by BAT, followed by a plateau lasting for a few 
thousands seconds, ending at the time Ta (following Willingale 
et al. 2007). This trend, named "Steep-Flat-Steep" (Tagliaferri et 
al 2005; Nousek et al 2006) has been interpreted in several ways 
(for a recent review see e.g. Zhang 2007). Furthermore, in nearly 
half of the bursts, X-ray flares, of relatively short duration At, i.e. 
At/t ~ 0.1, (e.g. Chincarini et al. 2007) are observed even several 
hours after the trigger. Considering X-ray flares in different GRBs, 
Lazzati, Perna & Begelman (2008) have shown that their average 
luminosity decays as t~ 5 ' 3 , similarly to what predicted following 
the mass accretion rate of fall-back material (see Chevalier 1989; 
Zhang, Woosley & Heger 2008). 

The optical light curves are also complex, but rarely track the 
X-ray flux behaviour (see e.g. Panaitescu et al. 2006, Panaitescu 
2007a, Panaitescu 2007b), suggesting a possible different origin. 
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For 10-15 per cent of bursts, precursor emission is detected, 
preceding the main event in some cases by hundreds of seconds. 
The energy contained in the precursors is comparable to that in the 
main event, and the spectra in the two phases are indistinguishable 
(Burlon et al. 2008), suggesting that they are produced by the same 
mechanism. 

Much theoretical effort has been made to understand the 
"Steep-Flat-Steep" behaviour, shown especially by the X-ray light 
curve. The initial steep decay is interpreted as "curvature" (or "high 
latitude") emission of the fireball: when the prompt ceases, the 
emission is dominated by radiation produced from parts of the fire- 
ball not exactly pointing at us (Zhang et al. 2006; Nousek et al. 
2006, but see e.g. Peer, Meszaros & Rees 2006 for an alternative 
interpretation). At later times, the relatively steep decay observed 
after Ta is generally explained as the "standard" forward (exter- 
nal) shock emission, namely as corresponding to the X-ray after- 
glow phase typically observed by BeppoSAX several hours after 
the burst trigger. There are however some alternative interpretations 
(see below). The most puzzling phase is the shallow, or plateau, 
one. Several models have been proposed, aiming at accounting not 
only for the shallow flux decay but also why it steepens at the break 
time Ta. The proposed alternative interpretations include: 
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• Energy injection. Zhang et al. (2006) propose that the shallow 
decay can be produced by a continuous, long-lasting, energy injec- 
tion into the forward external shock. There are at least two possibil- 
ities, depending on whether the central engine is long or short lived. 
A long-lived engine could have luminosity that smoothly decreases 
as L(t) oc t~ q with q < 1 (Zhang et al. 2006; Zhang & Meszaros 
2001), with Ta corresponding to the end of the energy injection 
phase. 

A short-lived engine (i.e. of duration comparable to that of the 
prompt phase) can produce shells with a steep power-law distribu- 
tion of r factors. Ta is determined by a cutoff in the Lorentz factor 
distribution. The two alternatives cannot be currently distinguished 
observationally. Both interpret the plateau as afterglow emission 
from a continuously refreshed shock. The required energetic (in 
bulk motion) largely exceeds what required to produce the prompt 
emission. 

• Reverse shocks. The shallow decay could be produced as syn- 
chrotron emission from the reverse external shock if the micro- 
physical parameters e c and eb are much larger than those in the for- 
ward shock. In this situation the ratio of the X-ray flux produced by 
the reverse and forward shocks would be dominated by the former. 
Along these lines Uhm & Beloborodov (2007) and Genet, Daigne 
& Mochkovitch (2007) suggested that the X-ray plateau emission 
is due to the reverse shock running into ejecta of relatively small 
(and decreasing) Lorentz factors. This requires an appropriate V 
distribution of the ejecta, besides the suppression of the X-ray flux 
produced by the forward shock. 

• Time dependent micro-physical parameters. If the relativistic 
electron distribution has a typical slope p ~ 2, the X-ray lumi- 
nosity is proportional to the bolometric luminosity (Lx = e c iboi). 
Since L bo i oc i -1 , a time evolution of e c oc t 1 ^ 2 would produce 
Lx oc t~ 1//2 , close to the observed plateau slopes. Ta is identified 
with the time e c reaches its maximum value (~ 0.1) (Ioka et al. 
2006). 

• Precursor fireball (Ioka et al. 2006). In this scenario a precur- 
sor, occurring 10 3 — 10 6 s before the main burst, generates a first 
fireball with low T, whose afterglow is too faint to be detected. The 
main burst then generates another, more powerful, fireball with a 
larger F. This second fireball, interacting with the first one, pro- 
duces the plateau phase. When the two fireballs merge and interact 
with the circum-burst medium, the standard afterglow sets in. 

• Up-scattering of forward shock photons. Panaitescu (2008) 
suggested that a relativistic shell would scatter protons produced 
by a forward shock located ahead of it. While this occurs also if the 
relativistic shell does not dissipate (bulk Compton), the process is 
more effective if dissipation (through, e.g. internal shocks) occurs, 
heating the electrons of the shell. The up-scattered component is 
expected to be more relevant in X-rays than in optical, thus over- 
shining the standard afterglow (i.e. forward shock) more easily in 
the X — ray band. 

• Geometrical models. If two co-aligned jets, with different 
opening angles, are observed at an angle 6 V within the wide cone, 
but outside the narrow one, the emission from the narrow jet would 
be visible only once it has decelerated. The observed light curve of 
the afterglow of the narrow jet would be flat before Ta, mimicking 
the observed plateau (see e.g. Racusin et al. 2008 for the case of 
GRB 0803 19B). In this model the time Ta is the time at which the 
Lorentz factor T of the narrow jet decreases to ~ 1/6 V . A some- 
what similar model is the off-beam model by Eichler & Granot 
(2006) in which the shallow phase represents the smooth peak of 
an afterglow observed off-axis. Here too Ta is identified by the 
time when the whole jet emission becomes visible. 



In the patchy shell model, Toma et al. (2006) propose that the 
early X-ray afterglow could be produced by an inhomogeneous 
jet of aperture angle ~0.1 rad composed by multiple sub-jets sub- 
tending a smaller aperture (# s _j < 0.01 rad). The latter ones are 
observable when their L decelerate to ~ 1/9 V . A shallow phase is 
ascribed to the superposition of the single sub-jet emissions, seen 
by an observer not exactly on-axis with any of them. When they 
merge due to sideways expansion the normal afterglow begins to 
dominate: this corresponds to the time Ta. Therefore the duration 
of the shallow phase depends on the (still uncertain) sideways ex- 
pansion velocity. 

• Dust scattering. Shao & Dai (2007) interpret the plateau as 
prompt X-ray flux scattered by dust grains located in the burst sur- 
roundings (within ~100 pc). This model has been recently ques- 
tioned by Shao et al. (2008) as it predicts a strong spectral soft- 
ening during the shallow decay phase, inconsistent with the data. 
Moreover, the large amount of dust required would imply an opti- 
cal extinction in excess of what observed. 

• Cannonballs. Dado, Dar & De Rujula (2005) propose that the 
entire steep-flat-steep behaviour of the X-ray light curve is due 
the sum of thermal bremsstrahlung and synchrotron emission from 
a cannonbal decelerating in the circum-burst medium. The time Ta 
would correspond to the start of the deceleration phase. 

• Ruffini et al. (2008) explain the early prompt and the 
steep-flat-steep phases within a unique scenario, a baryonic 
shell (fireshell) interacting with a non-homogeneous circum-burst 
medium. The emission process is thermal at all times. Once the 
fireshell reaches a region of low density, ~0.03-2 pc away, it would 
decelerate very slowly, giving origin to the plateau phase. 

It should be noted that the shallow phase is not ubiquitous: 
there are X-ray afterglows light curves where it is not detected. 
Therefore, any viable model should explain also the variety of the 
X-ray flux time behaviour. More importantly, all the above models 
propose that the X-ray shallow phase is due to afterglow emission 
(with the exception of the upscattering model by Panaitescu 2008). 
Thus, the same forward (or reverse) shock should produce optical 
radiation, which presumably would track the X-ray flux trend, in- 
cluding the shallow decay phase. This is not observed for some 
bursts which, therefore, challenge the above interpretation. 

Ghisellini et al. (2007) instead suggested that the plateau 
phase of the X-ray (and sometimes optical) emission corresponds 
to a "late prompt", namely due to the prolonged activity of the cen- 
tral engine (see also Lazzati & Perna 2007): after an early "stan- 
dard" prompt, the engine keeps producing shells of progressively 
lower power and bulk Lorentz factor for a long time (i.e. days). 
The dissipation during this and the early phases occurs at similar 
distances (close to the transparency radius). The reason for the shal- 
low decay phase, and for the break ending it, is that the F-factors of 
the late shells are monotonically decreasing, allowing to see an in- 
creasing portion of the emitting surface, until all of it is visible. The 
break at Ta occurs when F(t) — l/6j. In our scenario two inde- 
pendent emission components compete: the prevailing of the "late 
prompt" vs a standard afterglow emission at different times can ac- 
count for the variety of behaviours of X-ray and optical fluxes. 

In this work we thus try to model simultaneously both the X- 
ray and optical light curves as the sum of two components. The 
first one is the emission produced by the forward shock, accord- 
ing to the standard afterglow modelling. The second one is simply 
parametrised, spectrally and temporally. Though we refer to it as 
"late prompt" emission (which reflects our proposal), such a com- 
ponent could correspond to other interpretations. One of the aim 
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of our investigation is to find the constraints that a more physical 
model must satisfy to give origin to this "late prompt" component. 

Throughout this paper, a ACDM cosmology with Qm =0.3 
and SI a = h = 0.7 is adopted. 



2 THE SAMPLE 

All Swift bursts with known redshift, optical and X-ray follow up, 
as of end of March 2008, were considered. Among them, we se- 
lected GRBs for which an estimate of the optical extinction at the 
host site appeared in literature 0. This criterion is dictated by the 
need to determine reliable optical and X-ray intrinsic luminosities, 
in order to model their time dependent behaviour. The correspond- 
ing sample comprises 33 bursts. 

Information concerning these 33 GRBs are listed in Table 1, 
where we report: redshift, A^ Bt , optical spectral indices (3 (cor- 
rected for extinction), X-ray spectral indices (3x (again accounting 
for absorption), hydrogen column density N^ OBt (at the host) as 
determined by fitting the X-ray spectrum. 

It should be noted that usually A\? Bt is determined by re- 
quiring that the intrinsic optical spectrum is a power-law, and cor- 
recting the observed spectral curvature according to an extinction 
curve (the Small Magellanic Cloud one in most cases). Sometimes 
however, the requirement adopted was that the optical and X-ray 
data lie on the same functional curve (being it a single or a bro- 
ken power-law). When multiple choices were available, estimates 
based on the optical data alone were preferred as the X-ray flux 
could belong to a different spectral component. For details on the 
host frame dust absorption determination for each GRB see the ref- 
erences in Table 1 . 

In order to compare the behaviour of different bursts we de- 
reddened the observed optical fluxes taking into account both the 
GRB host dust and the Galactic (Schlegel, Finkbeiner & Davis 
1998) absorption along the line of sight. The reddening corrected 
fluxes have been then K-corrected and converted into monochro- 
matic luminosities through: 



L(v ) 



F{u ) 



(1) 



where vo is the central frequency of the photometric filter, g(l is the 
luminosity distance and (3 is the unabsorbed optical spectral index 
(see Table 1). 

The X-ray light curves were taken from the UK Swift Sci- 
ence Data Centrqj (see Evans et al. 2007 describing how the data 
were reduced). Also the X-ray 0.3-10 keV XRT light curves have 
been corrected for the combined effects of both host frame Ah and 
Galactic column densities, using the unabsorbed spectral index /3x 
obtained from the X-ray spectral analysis (see Table 1). The unab- 
sorbed 0.3-10 keV observer frame fluxes Fx have been converted 
to host frame 0.3-10 keV luminosities Lx as: 



£x = 



fx. 



(2) 



(1 + z) 1 ~^ 

For simplicity, we use the same /3x for the entire X-ray light curve, 
neglecting the sudden changes of /3x sometimes seen during X-ray 
flares, since the interpretation of the individual flares is beyond the 
aim of this work. The analysis has been carried on in the GRB host 



1 For one of them, GRB 070802, the photometric data set is not yet avail- 
able 

2 http://www.swift.ac.uk/xrt_curves/ 



time frame. We therefore rescale all the observed time intervals by 



3 THE MODEL 

As mentioned we assume that at all times the flux is the sum of two 
components: the first one due to synchrotron radiation produced by 
the standard forward shock caused by the fireball running into the 
circum-burst material; the second one is treated phenomenologi- 
cally, since its form/origin is not currently known, though it can 
be possibly ascribed to the extension in time of the early prompt 
emission. 



3.1 Forward shock component 

Following the analytical prescriptions of Panaitescu & Kumar 
(2000) the forward shock emission depends on the following pa- 
rameters: 

(i) Eo — the (isotropic equivalent) kinetic energy of the fireball 
after it has produced the early prompt radiation; 

(ii) To — the initial fireball bulk Lorentz factor. It controls the 
onset of the afterglow, but it does not influence the rest of the light 
curve. It is then rather undetermined when very early data are not 
available; 

(iii) no or M w /w w — no is the value of the circum-burst 
medium density if homogeneous, while Af w /i> w (wind mass loss 
rate over the wind velocity) determines the normalisation of the 
density in the wind case (oc R~ 2 ) profile; 

(iv) e — the "equipartition" parameter setting the fraction of 
the available energy responsible for electron acceleration; 

(v) 6b — the "equipartition" parameter parametrizing the frac- 
tion of the available energy which amplifies the magnetic field; 

(vi) p — the slope of the relativistic electron energy distribution, 
as injected at the shock. 

For simplicity, we assume that higher frequency of the afterglow 
synchrotron emission is beyond the X-ray range. These are 6 free 
parameters, if we consider no or M w /v w as a single one: in real- 
ity, the assumed homogeneous vs wind-like density profile can be 
considered as an additional degree of freedom. 

3.2 Late prompt component 

In the absence of a clear understanding of its origin, this component 
is parametrised with the only criterion of minimising the number 
of free parameters. This can be considered as a first step towards a 
more physical modelling of this second component. A subsequent 
analysis of the parameters distribution could help us in constraining 
possible theoretical ideas. A first attempt in this direction will be 
discussed in Sj5] 

The spectral shape - assumed to be constant in time - is de- 
scribed by a broken power-law: 

L L (v,t) = L (i)f _/3x ; 



L h (u,t) 



V > Vb 
v ^ i/ b , 



(3) 



where Lo is a normalisation constant. Lo is not treated as a free 
parameter by taking it as the 0.3-10 keV luminosity Llx of the 
late prompt emission at the time Ta: 



£-lx(7a 



Jo:. 



L h {v,T A )dv 



(4) 
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1.314 


0.72±0.06 


0.62±0.03 


0.87±0.08 


9.2±1.2 


52.85 


81.3 


Blo06, Pag07, Pag07 


061126 


1.1588 





0.93±0.02 


1.00±0.07 


11 ±0.7 


53.08 


70.8 


Per08a, Per08a, Per08a 


070110 


2.352 


0.08 


l.lzbO.l 


l.lzbO.l 


2.6zbl.l 


52.37 


88.4 


Jau07, Tro07, Tro07 


070125 


1.547 


0.11 ±0.4 


0.58±0.1 


l.lzbO.l 


2±1 


54.08 


65 


Fox07, Kan08, Upd08 


071003 


1.1 


0.209±0.08 


0.93±0.04 


1.14zb0.12 


1.1 ±0.4 


52.28 


150 


Per07, Per08b Per08b 


071010A 


0.98 


0.615±0.15 


0.76±0.25 


1.46±0.2 


17.4±4.5 


50.7 


6 


Pro07, Cov08a, Cov08a 


080310 


2.42 


0.1 ±0.05 


0.6 


0.9±0.2 


7.0±1 


52.49 


365 


Pro08, PeB08, Bea08 


0803 19B 


0.937 


0.07±0.06 


0.33±0.04 


0.814=b0.013 


1.87±0.13 


53.27 


50 


Vre08, BI0O8, BI0O8 



Table 1. The sample. For all bursts we report information taken from the literature (see the references), namely: redshift, optical extinction and hydrogen 
column density at the host (/ly ost andN^ st , respectively), and the optical and X-ray indices found after de-absorbing. E^ so is in the 15-150 keV band, 
not K-corrected. T90 is in seconds, from the Swift catalogue (http : //swift ■ gsf c ■ nasa ■ gov/docs/swift/archive/grb_table ■ html 7) . 
References: Ber05a: Berger et al. (2005a); Sti05: Still et al. (2005); Per05: Perri et al. (2005); Fyn05a: Fynbo et al. (2005a); Kan08: Kami et al. (2008); 
Cus06: Cusumano et al. (2006); Fyn06: Fynbo et al. (2005b); Wat06: Watson et al. (2006); Ber05b: Berger et al. (2005b); DUP07: importantde Ugarte Postigo 
(2007); Cap07: Capalbi et al. (2007); Cen05: Cenko et al. (2005); Hol07: Holland et al. (2007); Man07a: Mangano et al. (2007a); Fol05: Foley et al. (2005); 
Blu06: Blustin et al. (2006), Che05: Chen et al. (2005); Pan06: Pandey et al. (2006); Per06: Perri et al. (2007); DeP07: de Pasquale et al. (2007); Fyn05d: 
Fynbo et al. (2005d); Sch07: Schady et al. (2007); Oat07: Oates et al. (2007); Pro05: Prochaska et al. (2005); Cen06a: Cenko et al. (2006a); Pag05: Page et al. 
(2005); Fyn05f: Fynbo et al. (2008f); Sol07: Sollerman et al., (2007); Jak05a: Jakobsson et al. (2005a); Ken05: Kennea et al. (2005); Hil05: Hill et al. (2005); 
Gui07: Guidorzi et al. (2007); Cen06b: Cenko et al. (2006b); Rom06: Romano et al. (2006); Fyn06: Fynbo et al. (2006a); Mor06: Morris et al. (2006); Cuc06: 
Cucchiara Fox & Berger (2006); Cur07b: Curran et al. (2007b); Pro06: Prochaska et al. (2006); E1106: Ellison et al. (2006); Fal06: Falcone et al. (2006); 
BI0O6: Bloom et al. (2006); Ber06: Berger & Gladders (2006); Tho08: Thone et al. (2008), Cam06: Campana et al. (2006a); Pri06: Price Berger & Fox (2006); 
Man07b: Mangano et al., (2007b); Th606: Thone et al., (2006); Gru07: Grupe et al. (2007); Fug06: Fugazza et al. (2006); Gru06: Grupe et al. (2006); R0IO6: 
Rol et al. (2006); Eva06: Evans et al., (2006); Fyn06c: Fynbo et al. (2006c); RuV07: Ruiz-Velasco et al., (2007); Osi06: Osip Chen & Prochaska (2006); 
Mun07: Mundell et al. (2007); BI0O6: Bloom Perley & Chen (2006), Pag07: Page et al. (2007); Per08a: Perley et al. (2008a); Jau07: Jaunsen et al. (2007); 
Tro07: Troja et al. (2007); Fox et al. (2007); Upd08: Updike et al. (2008); Per07: Perley et al. (2007); Per08b: Perley et al. (2008b); Pro07: Prochaska et al. 
(2007); Cov08: Covino et al. (2008a); Pro08: Bea08: Beardmore et al. (2008), PeB08: Perley & Bloom (2008a); Prochaska et al. (2008); Vre08: Vreeswijk et 
al. (2008); BI0O8: Bloom et al. (2008). 



with v in keV. Again for simplicity we assume that any cut-off 
frequency, at high as well as at low energies, is outside the IR- 
optical/X — ray frequency range. 

The temporal parameters, described by the flat and steep decay 
indices, aa and a st respectively, and the time Ta at which the two 
behaviours join, are assumed to be described by a smooth broken 
power-law: 



L L (>,t) = L h (v,T A ) 



(*/*a) 



1 + {t/t A ) a ^- aci ' 



(5) 



To summarise, the free parameters reproducing the late prompt 
emission are: 



(i) /3x — the spectral index of the late prompt emission in X- 
rays; 

(ii) P — the spectral index of the late prompt emission in the 
IR-optical; 

(iii) fb — the break frequency between the optical and the X- 
rays; 
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(iv) Llx(Ta) — the 0.3-10 keV luminosity of the late prompt 
emission at the time Ta; 

(v) afl — the decay index for the shallow phase, before Ta 

(vi) a s t — the decay index for the steep phase, after Ta; 

(vii) Ta — the time when the shallow phase ends. 

These are 7 free parameters. It is worth stressing that, despite of 
their number, these are rather well constrained by observations. 
When the late prompt emission dominates, a&, a st , Ta can be di- 
rectly determined as well as one spectral index (usually /3x, since 
the late prompt emission is usually dominating in the X-ray range). 
Some degeneracy is present between V\> and /3 , both of which con- 
trol the importance of the optical flux due to the late prompt compo- 
nent: the same optical flux can for instance be reproduced assuming 
a steeper (flatter) (3 and a larger (smaller) v^, as the ratio between 
the 0.3-10 keV X-ray luminosity and the v L(v ) optical luminos- 
ity of the late prompt is proportional to u^ x ~ ■ 

3.3 Caveats 

As our treatment is necessarily simplified, simply parametrising 
the late prompt emission, we analyse below the most important (or 
drastic) assumptions, trying to outline their effects. 

• The afterglow calculations are based on the prescriptions by 
Panaitescu & Kumar (2000). In their analytical treatment the cur- 
vature of the emitting shell is neglected. The inclusion of the time 
delay between the emission times of photons received at any ob- 
server time would smooth out any relatively sharp feature of the 
light curve (especially when the injection or cooling frequency 
crosses the considered band). However the derived light curves are 
sufficiently accurate for the purposes of the present work. 

• Almost all of the calculations of the afterglow light curves as- 
sume that e c and eb are constant in time. This is likely to be just 
a rough approximation, since the physical conditions at the shock 
front change in time (r as well as the density measured in the co- 
moving frame do change). As such a temporal dependence is not 
known or predicted, we are forced to adopt this simplification. 

• The afterglow emission is assumed to be isotropic, therefore 
no jet breaks can be reproduced in the calculated light curves. 

• The spectrum of the late prompt emission is assumed to be 
constant in time, in the observer frame. This is likely to be the most 
critical approximation, adopted just to minimise the number of free 
parameters. One might speculate that if this component originate 
by shells with decreasing bulk Lorentz factor (as in the models by 
Uhm & Beloborodov 2007, Genet, Daigne & Mochkovitch 2007 
and Ghisellini et al. 2007), then it is likely that the observed break 
frequency v\, would also decrease in time (if constant in the co- 
moving frame). While this would not affect the X-ray light curves 
(if i/b is below the X-ray window even at early times), the optical 
emission would become relatively more important as time goes on. 
For instance, a plateau in the X-rays could correspond to a rising 
optical light curve. This suggests a possible observational test. As- 
sume to select a burst in which both the optical and the X-ray light 
curves are dominated by the late prompt emission. If decreases 
in time, we should see two effects. First, the optical plateau should 
be shallower than the X-ray one (since the X-ray to optical flux 
ratio decreases as v^~^°). Secondly, when v\> crosses the optical 
band, we should see a spectral steepening, since the decreasing 
acts as a cooling break. After v\, has crossed the optical band, the 
optical and the X-ray fluxes should lie on the same power-law. 

• The low and high frequency cut-offs of both the afterglow and 
late prompt emission have been neglected as free parameters. The 



late emission spectrum might have a high frequency cut-off in the 
X-ray band. Given the current status of the X-ray observations, 
that do not detect such a cut-off, this simplification is reasonable. 

• The late prompt emission is assumed to last forever, while, of 
course, it will die away after some time. This may happen, however, 
at very late times, when any X-ray or optical observations are not 
any longer feasible or when the GRB emission cannot be detectable 
(in the optical, emission can be dominated by the host galaxy or 
sometimes by a supernova associated to the burst). 

• Flares, re-brightenings and/or bumps in the light curve are 
not acconted for. In our scenario, these are separated components, 
though in practice, their presence makes the choice of what data 
points to "fit" a bit subjective. 

A final remark. Due to the above caveats, the values of the 
parameters for a single source may be subject to rather large uncer- 
tainties. In this sense the distributions of parameter values are much 
more meaningful. We could badly model an individual source, but 
the general conclusions could be right, if some coherence is found 
for the parameters of the entire sample. 



4 RESULTS 

Figs. [JJ{9] show the X-ray and optical light curves of the 33 GRBs 
together with the results of the modelling: dotted lines refer to the 
late prompt emission, dashed ones to the afterglow component and 
the solid lines to their sum. 

The parameters inferred from the modelling of the light curves 
are reported in Table(2] together with a tentative classification of the 
bursts according to the dominant contribution: "A" stands for after- 
glow, "L" for late prompt, and "X" and "O" refer for X-ray and 
optical, respectively. For instance, XL-OA indicate that the X-ray 
flux is dominated by the late prompt, and the optical by the after- 
glow. When both type of emissions are comparable we use "M", 
for mix. This also comprises the case when one component dom- 
inates in one time interval, and the other in another time interval. 
The number of bursts which can be described within these cate- 
gories is summarised in Table [4] The X-ray flux is dominated by 
the late prompt emission or a mixture of late prompt and afterglow 
for the majority of GRBs, the opposite being true for the optical 
emission. Out of our 33 events, the most common cases are XM- 
OA (10 GRBs, namely a mix in the X-rays and afterglow in the 
optical) and XL-OM (8 GRBs, namely late prompt in the X-rays 
and a mix in the optical). 

The overall result is that both components have comparable 
relevance in most cases. This is can be seen as a direct consequence 
of the different slopes of the light curves: since the late prompt 
emission is flatter than the afterglow up to Ta, and often steeper 
after this time, it is likely that the late prompt emission dominates 
or contributes around Ta even in the optical. Conversely, the af- 
terglow may dominate or be important at very early and very late 
times (if there is no jet break). In other words the similar contribu- 
tion of both components is the cause of the complex X-ray-optical 
behaviour observed. 



4.1 Distribution of parameters 

Figs. [TO] and QT] show the distribution of all our input parameters. 
For comparison, in these figures we report the values found by 
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Table 2. Input parameters for the afterglow component (columns 2-7) and for the late prompt emission (columns 8-14). Col 1 : Burst Id; Col 2: Fireball kinetic 
energy (after the early prompt emission, in units of 10 53 erg); Col 3: Initial bulk Lorenz factor; Col 4: density of circum-burst medium: values equal or larger 
than 1 are for a homogeneous density; values much smaller than 1 correspond to a wind like profile; the listed value is M w /v w , where M w is the mass loss 
rate in A/Q/yr and v w is the wind velocity in km s^ 1 . Col. 5 and 6: equipartition parameters e c and e^; Col 7: slope of the assumed relativistic electron 
distribution; Col. 8: spectral break of the late prompt emission (in Hz); Col. 9 and 10: high and low energy spectral indices of the late prompt emission; Col. 
1 1 and 12: decay slopes of the late prompt emission, before and after Ta listed in Col. 13 (in sec); Col. 14: luminosity (in units of 10 45 erg s^ 1 ) in the 0.3-10 
keV energy range of the late prompt emission, at the time Ta; Col. 15: Burst classification (see text). 



Panaitescu & Kumar (2002) for 10 pre-Swift burst£j| This is to be 
expected, since in our interpretation the X-ray luminosity in the 
majority of cases is not produced by the afterglow, which is thus 
less energetic. 

Most (25 out of 33) afterglows can be consistently described 
by the interaction of the fireball with a homogeneous medium. This 

3 Panaitescu & Kumar (2002) give the collimation corrected value for the 
isotropic kinetic energy of the fireball after the early prompt phase. We have 
then divided this value by (1 — cosfij) to get the isotropically equivalent 
value of Eq to be compared with the values found for our bursts. Note 
that To does not affect the properties of the afterglow after its onset, and is 
therefore not an important parameter for Panaitescu & Kumar (2002), who 
are fitting data taken much later than the afterglow onset (with the exception 
of GRB 990123). The afterglow parameters found for our bursts are rather 
standard, being similar to the ones obtained by Panaitescu & Kumar (2002) 
(see also Panaitescu & Kumar 2001a and 2001b). The distribution of the 
circum-burst density no is narrower for the bursts in our sample, while the 
distributions of e c and eg are centred on smaller values. 



is especially the case when the optical light curve indicate the onset 
of the afterglow itself (i.e. a very early rising phase), that cannot be 
reproduced with a wind-like density profile. The latter in fact pro- 
duces almost flat optical light curves in the early phases. The ho- 
mogeneous densities are very narrowly distributed around a mean 
value of (no) ~ 3 cm -3 . For 8 GRBs (see Tab. O a better mod- 
elling can be achieved invoking a wind-like density profile. All but 
one of these 8 bursts can be modelled with a value of the ratio of 
the mass loss rate and the wind velocity of M w /« w = 10 -8 (Mq 
yr _1 )/(km s -1 ) that can correspond to M w = 10 -5 Mq yr _1 and 
b w = 10 3 km s _1 . The remaining burst require half of this value. 
Similarly to what had been found by Panaitescu & Kumar (2002) 
the afterglow parameters distribution are quite broad, i.e. they do 
not cluster around typical values. Exceptions are the density no and 
the bulk Lorentz factor To . 

Also the distributions of some late prompt parameters (i.e. Ta, 
Lt a and i^b) are rather broad, while (3 and the temporal slopes qa 
and a st are more narrowly distributed. The values of Ta range from 
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Figure 1. Figs 1-9. X-ray (in grey) and optical (different symbols, as la- 
belled) light curves. Lines indicate the model fitting: afterglow (dashed), 
late prompt (dotted) and their sum (solid). Black lines refer to the X-rays, 
light grey (red in the electronic version) for the optical. The vertical line 
(and shaded band) correspond to the rest frame jet break times (and their 
3<t uncertainty). Grey lines and stripes correspond to jet break times as re- 
ported in the literature (references are listed in Ghirlanda et al. 2007), light 
grey (yellow in the electronic version) lines and stripes refer to jet times 
expected if the burst followed the Ghirlanda relation. These are shown only 
for bursts with measured -E pea k> me peak energy of the prompt emission. 
References for the optical data can be found in the appendix. 
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Figure 2. Same as in Fig. 1. 
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Figure 5. Same as in Fig. 1. 



Figure 6. Same as in Fig. 1. 
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Figure 7. Same as in Fig. 1. Figure 8. Same as in Fig. 1. 
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Table 3. Number of sources dominated by different components: XA (OA): 
X-ray and optical flux dominated by the Afterglow emission; XL (OL): 
X-ray and optical flux dominated by the late prompt emission; XM (OM): 
X-ray and optical fluxes where the late prompt and afterglow emission are 
relevant. 

10 2 to 10 4 s or more (in the rest frame), and are (anti-)correlated 
with the late prompt luminosity at Ta, as shown in Fig. [12] This 
confirm the correlation found by Dainotti, Cardone & Capozziello 
(2008). This results in a narrow distribution of T\Lt a (Fig. II It. 

The distributions of j3 a and v a must be taken with caution, 
since the model fixes only their combination, and only in a few 
GRBs they can be constrained separately (i.e. when the optical light 
curve is dominated by the late prompt emission and the spectral 
index during this phase is known). 

The distribution of a st is intriguing, since it is centred around 
a mean value of 1.6. This is very close to 5/3, the predicted decay 
of the accretion rate of fall-back material (see also §5 where this 
point is discussed in more depth). The values of aa cluster around 
0. 

In Fig.QT]we show the distribution of TaLt a , and in Fig. 1 131 
we show TaLt a as a function of Ei so . The two quantities are cor- 
related (albeit poorly) and the energy contained in the late prompt 
emission (of which Ta£t a is a proxy) is at most comparable with 
Ei so . More frequently TaLt a is one or two orders of magnitude 
smaller than Ei so , in agreement with the findings by Willingale et 
al. (2007). 

Fig. QT] shows also the distribution of Ei so /[Eo + -Eiso]- This 
ratio represents 77, the fraction of the total energy of the fireball re- 
quired to produce the observed early prompt radiation. In Fig. [14] 
this fraction is shown as a function of -Ei so - Although there is a 
weak positive correlation, the mean value is well defined and cor- 
responds to 77 ~ 0.1. 

4.2 Jet breaks 

A currently hot debate concerns the absence of jet breaks in the 
light curves of GRB afterglows. In the scenario we propose the light 
curve comprises two components of which only the afterglow one 
should present a jet break (at tj ). It follows that jet breaks should be 
more often detectable in the optical, rather than being achromatic, 
and the after-break slopes may be shallower than predicted by the 
closure relations. 

No jet breaks — When the flux is dominated by the late prompt 
emission in both the optical and the X-ray bands, jet breaks may 



Time [s rest frame] 




Time [s rest frame] 
Figure 9. same as in Fig. 1 . 

become unobservable. The late prompt emission (at least after a few 
thousand seconds) does so for 6 GRBs of the sample (namely GRB 
050319, GRB 050408, GRB 060614, GRB 060729, GRB 061126 
and GRB 071003). Therefore, for these bursts, no jet break is pre- 
dicted to be visible if the late prompt light curve continues unbro- 
ken for a long time - if the late prompt component instead breaks, 
we might erroneously interpret this as a jet break. 
Achromatic jet breaks — Viceversa, an achromatic jet break 
should be observed when both in the optical and X-ray light curves 
the afterglow emission prevails, at least when the jet break is likely 
to occur. 16 GRBs of the sample could show such an achromatic 
break (GRB 050318, GRB 050401, GRB 0504 16A, GRB 050802, 
GRB 050820A, GRB 050824, GRB 060512, GRB 060904B, GRB 
060908, GRB 060927, GRB 061121, GRB 070110, GRB 070125, 
GRB 071010A, GRB 080310, GRB 080319B). Emission in several 
of these bursts - although dominated by afterglow emission, espe- 
cially in the X-ray band, at late times - still comprises a relevant 
contribution from the late prompt component. Therefore the steep- 
ening of their light curve after tj should be shallower than what the 
standard afterglow theory predicts. 

Chromatic jet breaks — When the late prompt is dominating in 
one band, and the afterglow in the other, a jet break should be vis- 
ible only in the afterglow-dominated band. According to our find- 
ings a jet break could be present in the optical but not in the X- 
rays band in 9 GRBs (GRB 050525 A, GRB 050730, GRB 050801, 
GRB 050922C, GRB 060124, GRB 060206, GRB 060418, GRB 
060526 and GRB 061007). Instead, 2 GRBs (GRB 051111 and 
GRB 060210), could show a jet break in X-rays but not in the op- 
tical. 
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Figure 10. Top 4 panels: distribution of the values of the micro-physical 
parameters e c , eg, homogeneous density no, and electron slope p. Bottom 
4 panels: distribution of the isotropically equivalent initial kinetic energy 
Eq, bulk Lorentz factor To, break frequency and optical spectral in- 
dex for the late prompt emission /3 . The hatched areas correspond to the 
distribution of parameters found by Panaitescu & Kumar (2002) fitting the 
afterglow of 10 pre-Swift bursts. They are shown for comparison. 
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Figure 11. Top 4 panels: distributions of the decay indices of the late prompt 
emission, and a s t, of Ta and of the 0.3-10 keV luminosity at the time 
T\. Bottom 4 panels: distributions of the isotropic energy i? 7i i so of the 
early prompt radiation, of the ratio E-, ; so / (E~ t i so + Eq), which provides 
an estimate of the efficiency of the prompt emission; of the energy T A Lj- A , 
and of the ratio T\Lt a / '-Esso- 



in Figs.[T}[9]we indicate the time at which a jet break has been 
reported to be detected or the time at which a jet break is expected 
to be seen if the burst were to follow the -E pca k — (Ghirlanda) re- 
lation (Ghirlanda, Ghisellini & Lazzati 2004, updated in Ghirlanda 
et al. 2007) (see the figure caption). The latter ones are estimated 
only for bursts with measured -E pca k, the peak energy of the vF v 
spectrum of the proper prompt emission. We found no contradic- 
tory cases (i.e. an observed jet break occurring in a late prompt- 
dominated GRB), except for GRB 060614. 

There are some additional bursts for which the presence of a 
jet break has been claimed in the literature. For instance, in GRB 
050319 Cusumano et al. (2006) suggest that the break in the X- 
ray light curve at 27,000 s (observed time) could be a jet break, 
but also discuss the problems with this interpretation due to the 



unusual pre- and post-break slopes. In our scheme, the observed 
break simply corresponds to 7a- 

For GRB 050730, Pandey et al. (2006) consider the change 
of slopes at ~ 0.1 d (observed time) in the optical light curve as 
indicative of a jet break. In our interpretation, instead, the change 
of the flux decay slope is due to the late prompt emission providing 
a relevant contribution after ~ 3 x 10 3 (rest frame time). 

Malesani et al. (2007) claim the presence of a possible jet 
break in the optical light curve of GRB 070110, at ~ 5 days (ob- 
served time). According to our findings, this can indeed be a jet 
break that should also be visible in X-rays. 

In the light curves examined here, there are also a few exam- 
ples of slope changes that could be jet breaks, but for which we 
could not find any report in the literature. The optical light curve of 



© 0000 RAS, MNRAS 000, 000-000 



A unifying view of Gamma-Ray Burst Afterglows 13 



10 51 

10 50 

10 49 |- 
10 4B %■ 
10 47 
10 46 t- 

10 45 r 



< 



10 52 
10 61 
10 50 k 



10 54 r 
10 53 |- 
10 52 
10 51 
10 50 L- 



o 

°°c 
o c 



oo 



I I x 



o 

I I 1 1 III 



8° 

o 



o 
oo 



I I I T III i 



O o 



o 

OO^O o 



o 
oo 



oo 



o 



o o 
— I — I I 1 1 1 III 



oo 



§)0 ° oo' 

® ° °°oo° 
oo 



,Qo 



10 100 1000 10 4 

T A [s, rest frame] 



io 5 



Figure 12. The luminosity of the late prompt emission Lt a (in erg s ) 
at Ta, the corresponding energy TaLt a (in erg) and the isotropic energy 
^7, iso (in el 'g) as functions of T&. Note that Ly A anti-con'elates with Ta, 
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Figure 13. Energy of the late prompt emission, estimated as T\Lt a , as a 
function of the isotropic energy of the prompt emission, £L i so . The dashed 
line corresponds to the least square fit, [T\Lx A ] oc E®'f® (chance proba- 
bility P = 2 X 1CT 7 , excluding the outlier GRB 070125). 



GRB 060206 may be one of such cases (see the last optical point in 
Fig-EJ. For this GRB the presence of the jet break is expected only 
in the optical, since the X-rays are dominated by the late prompt 
component. Note that the corresponding tj would make this burst 
consistent with the Ghirlanda relation (see the vertical grey line in 
Fig-EJ- Another example is visible in the X-ray flux decay of GRB 
061121, at ~ 10 J s (rest frame, see Fig.0. Unfortunately, there 
are no optical data at this late time to confirm it. Again, if this is a 
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Figure 14. Top panel: the kinetic energy Eo after the prompt emis- 
sion as a function of EL j so . The dashed line is a least square fit, yield- 
ing Eo oc E®-f^ o (chance probability P ~ 10 -3 , excluding GRB 
070125). Bottom panel: the efficiency of the prompt emission estimated as 
E-y, iso/ (E-y^so+E q) as Si function of -E 7 ,i so . There seems to be weak cor- 
relation, in the sense that weaker bursts would have the smaller efficiency. 
See the corresponding distribution in Fig. II II Here and in the other figures, 
the plotted values of -E 7i i so are neither bolometric nor K-corrected, but 
refer to the observed 15-150 keV range. 



jet break, the burst would be consistent with the Ghirlanda relation 
(see the vertical grey line in Fig. [7}. Also in GRB 071010 there 
could be a jet break in optical, after ~ 10 5 s (rest frame, see Fig. 
[8} but its interpretation is difficult because of an optical/X-ray flare 
occurring just before. Finally, for GRB 080310, a steepening of the 
optical light curve after ~ 10 s (rest frame, see Fig. [9} could be a 
jet break, as also supported by a steepening also in the X-ray light 
curve, that is (marginally) dominated by the afterglow component. 

We plan to discuss in more detail these possible jet breaks in 
a forthcoming paper (Nardini et al., in preparation). 



4.3 Prompt and afterglow energetics 

As the X-ray luminosity ix is found to be often dominated by 
the late prompt emission, it does not provide a proxy for the after- 
glow bolometric luminosity. Since Lx exceeds what observed in the 
other spectral bands, the estimated luminosities and total energetics 
produced by the afterglow are radically smaller than what simply 
inferred from Lx- 

This exacerbates the problem of understanding why the early 
prompt emission is larger than the afterglow one, if the former is 
dissipated in internal shocks. In fact, while in external shocks, be- 
lieved to be responsible for the afterglow, a fraction of the whole 
fireball kinetic energy is available, in internal shocks only a fraction 
of the relative kinetic energy between two colliding shells can be 
dissipated as radiation. If such fractions are similar, the "bolomet- 
ric afterglow fluence" is expected to be a factor ~10 larger than the 
bolometric early prompt fluence. The opposite is observed, and the 
discrepancy is more extreme if Lx provides only an upper limit to 
the afterglow contribution, as in our interpretation. 
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Figure 15. The energetics of the afterglow component, estimated as e e Eo, 
as a function of: (top panel) -E 7i i so , the energetics of the prompt emission 
as measured in the 15-150 keV band (rest frame) (top panel) - the dashed 
line corresponds to equal values; (bottom panel) -Eiato, the energetics of the 
late prompt emission, as measured in the (rest frame) 0.3-10 keV band and 
approximated by T a Lt a ■ 



Bearing in mind that it is often dangerous to claim correlations 
between luminosities or energetics, since both quantities are func- 
tion of redshift, we can compare Fig.[l2]wifh Fig.[T3] It can be seen 
that the correlation between the late prompt energetics measured by 
TaLt a and i5 7 ,i so is stronger than the correlation between the ki- 
netic energy (after the early prompt) Eq and E lt i so . A least square 
fit yields [T a L Ta ] oc £°;f s 6 (chance probability P = 2 x 10~ 7 ), 
and E a oc £°;f s 2 (chance probability P ~ 10" 3 ). If the T a L Ta - 
E-y^so relation is not a mere product of the common redshift depen- 
dence (which however should also affect the Eo-E lt i so relation) 
this suggests that the early and the late prompt phases of emission 
are related. 

In Fig. [15] (top panel) e c Eo (which can be considered as an 
upper limit to the bolometric afterglow luminosity) is compared to 
i? 7 ,iao, the energetic of the prompt emission as measured in the 15- 



150 keV band (rest frame). _E 7i i so exceeds the afterglow energetics 
by almost two orders of magnitudes. In the bottom panel of the 
same figure e c Eo is plotted against the energetics of the late prompt 
emission £a a tc, approximated by the quantity TaZ/t a • These quan- 
tities do not correlate, suggesting that they are two separated com- 
ponents. 

To summarise: all indications gathered from the analysis of 
the energetics suggest that what we have called "late prompt emis- 
sion" is a phenomenon not related to the afterglow, but it is more 
connected to the same engine producing the early prompt. Further- 
more, the energetics associated to the afterglow emission is on av- 
erage a small fraction of the total energy of the burst. 



5 DISCUSSION AND CONCLUSIONS 

The proposed scheme appears to be suitable to account for the di- 
versity of the optical and X-ray light curves of GRBs, at the ex- 
pense of introducing, besides the standard afterglow emission re- 
sulting from the external (forward) shock, another component. This 
has been simply parametrised with 7 free parameters. The distri- 
butions of these parameters are not particularly clustered around 
mean values, except for the time decay slopes qa and a at (see 
below). However, this should not be taken as a potential problem 
for the proposed idea, since even the well established afterglow 
model, when applied to the optical and X-ray afterglows of pre- 
Swift GRBs, yield broad parameter distributions (see Fig. 1 101 and 
Panaitescu & Kumar 2002). 

Our phenomenological approach should be considered as a 
first step towards the construction of a convincing physical model. 
As discussed in the introduction, there has been already a bloom- 
ing of theoretical ideas, but a general consensus has not yet been 
reached. Our findings can shed some light and help to discriminate 
among the different proposals. 

As an illustrative example, the proposed scenario can be con- 
trasted with the alternative idea that GRBs are characterised by two 
jets with different opening angles (see the Introduction). In the lat- 
ter interpretation, if the line of sight lies within the wide cone but 
outside the narrow one, the emission from the narrow jet will be ob- 
servable when T has decreased to T ~ 1 /8 V and the corresponding 
afterglow light curve can reproduce the flat-steep-flat behaviour 
and present a break (at Ta). However, it is hard to explain why the 
flat-steep-flat trend is not observed in the optical, as in a (narrow 
jet) afterglow the optical and X-ray fluxes should temporally track 
each other. The "late prompt" scenario appears to provide a better 
interpretation of the data. 

Within the proposed scheme, some light can be shed on the 
puzzling issue about jet breaks. 

They can be achromatic if the afterglow component is ob- 
served in different bands. These may not be the case for several 
bursts. Furthermore, there are a few in which the late prompt emis- 
sion, dominating both the optical and X-ray flux, hides jet breaks at 
all times. For these bursts a break at the time Ta is visible in both 
the X-ray and optical bands, and it can be erroneously be taken 
as a jet break, since it is achromatic. Only a densely sampled light 
curve in both bands can help to discriminate the presence of a real 
jet break. Note that the bursts in our sample, selected for having 
an estimate of extinction in the host, are better sampled optically 
than the rest of the bursts, for which it would be difficult to reliably 
estimate the relative importance of the two components. 

Another relevant consequence of our scenario is that, even 
when a jet break is observed, the after-break light curve can be 
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Figure 16. The light curves of all the 33 GRBs in the X-rays (left panel) and optical (right panel). For comparison, the dashed lines correspond to t~ 5 / 4 and 
t~ 5 / 3 , as labelled. Especially in the X-rays, the luminosity profile seems to be flatter than t 5 / 3 and closer to a t~ 5 / 4 decay. However, this behaviour is due 
to the contribution in some GRBs of the afterglow emission at late times, flattening the overall light curve. See Fig. ll7l for comparison. 



flatter than predicted, since the late prompt flux can contribute af- 
ter tj . This implies that the so-called closure relations, linking the 
flux decay slopes before and after the break with the spectral in- 
dex, should be taken with care. In this respect, it is worth stressing 
that the closure relations, and in general the simplified afterglow 
scenario predicting them, treat the micro-physical parameters e c 
and £b as constant in time. We adopt the same simplification, but 
this might become a crucial issue once we will have a convincing 
physical interpretation predicting their time behaviour. 

^.Frorn our modelling the steep decay of the late prompt emis- 
sion can be described by a power-law with slope a st ~ 1.6. This 
is intriguingly similar to the time dependence of the mass accretion 
rate during the fall-back phase, and to the average decay of the 
X-ray flare luminosity, as analysed by Lazzati, Perna & Begelman 
(2008). This is not the average decay slope observed: the X-ray and 
optical light curves are flatter than L(t) oc i -5 / 3 (see Fig. 1161), but 
this is due to the contribution, especially at early and late times, 
and in the X-rays, of the afterglow contribution. Fig. [T7] shows 
the results of our light curve modelling for the optical and X-rays 
bands. The late prompt light curves are indeed steeper, on aver- 
age, than the sum of the two components that reproduce the data. 
We consider this as a main result of our analysis, because it sug- 
gests that the late prompt emission can be interpreted as due to the 
late time accretion onto a black hole of fall-back mass, namely 
material that failed to reach the escape velocity from the explod- 
ing progenitor star, and falls-back. According to analytical results 
(Chevalier 1989) and numerical simulations (e.g. Zhang, Woosley 
& Heger 2008), the accretion rate decreases in time as t -5 ' 3 , and 
can continue for weeks, enough to sustain late prompt emission 



even at very late times. Our finding also agrees with that obtained 
by Lazzati, Perna & Begelman (2008) by analysing X-ray flares. 
They found that the average luminosity of X-ray flares, for a sam- 
ple of GRBs with known redshift, also decays like t~ 5 ^ 3 . Such 
an agreement then suggests that both the X-ray flares and the late 
prompt emission have a common origin, related to the accretion of 
the fall-back material. It remains to be explained why this phase is 
observed after Ta that in some cases can be as long as 10 4 sec or 
more, while the simulations predict a quasi-constant accretion rate 
for 10 2 -10 3 s (MacFadyen, Woosley & Heger 2001). There are at 
least two possibilities. The first one is suggested by the simulations 
of Zhang, Woosley & Heger (2008) (see their Fig. 2) which include 
the effect of the reverse shock running through the fall-back mate- 
rial. When the reverse shock reaches the inner base the material is 
slowed down, and thus the accretion rate is enhanced. The asymp- 
totic i -5 / 3 phase can thus be delayed. The second possibility has 
been suggested by Ghisellini et al. (2007): even if the total flux 
produced by the late prompt phase is decaying at the rate t~ 5 ^ 3 , a 
decreasing T implies that the observed emission comes from an in- 
creasing surface (oc 1 /r 2 ), making the observed decay flatter than 
t~ 5/3 , until, at T A , r ~ l/0j. After T A the whole emitting surface 
contributes to the detected flux, and the flux decreases as i -5//3 . 

Fig. [T7] shows that the sum of the late prompt and afterglow 
emission makes the optical fluxes to cluster. This occurs because 
the late prompt emission - though usually not dominant in the opti- 
cal - narrows the distribution of the optical luminosities at a given 
time. The vertical dotted line in the figure corresponds to the time 
(12 hours)) at which Nardini et al. (2006, 2008), Liang & Zhang 
(2006), Kann, Klose & Zeh (2006), found a remarkable clustering 
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Figure 17. The light curves, as inferred from the modelling, for the 33 GRBs in the X-ray (top panels) and optical (bottom panels) bands. The late prompt 
(left panels), afterglow (middle panels) and total (left panels) emission are shown. The vertical dotted lines correspond to 12 hours. Note that the total optical 
luminosity at 12 hours is more clustered than the late prompt and afterglow luminosities. The dashed lines in the left panels correspond to L oc t~ 5 / 3 , while 
in the middle and right panels also decays L oc t~ 5 / 4 are shown for reference. 



of the optical luminosities around two well separated values. How- 
ever the total optical luminosities (right bottom panel of Fig. |17t 
are more dispersed than the afterglow ones (middle bottom panel). 

Our scenario makes some predictions and calls for some con- 
sistency checks. We are analysing in more details the data for the 
GRBs of the present sample to find confirmation of and/or prob- 
lems with our scenario, and the results will be presented in a forth- 
coming paper (Nardini et al., in preparation). The main obvious 
prediction concerns the presence or absence of jet breaks. When 
both the X-ray and optical light curves are late-prompt dominated, 
no jet break should be seen. Viceversa, when they are both dom- 
inated by the afterglow emission, an achromatic jet break is ex- 
pected (even if the after-break slope may be shallower). A chro- 
matic jet break should be observed when only one of the two spec- 
tral bands is dominated by the afterglow flux. We stress that well 
sampled data are required to reliably assess the relative contribution 
of the two components. 

A further general prediction concerns the spectral shape of the 
late prompt flux. In order for it to be negligible with respect to 
the the afterglow optical emission, there must be a spectral break 
between the optical and X-ray bands, and the slope below the break 
should be rather flat. From our modelling some GRBs require that 



either this slope is extremely hard or there is a break within the X- 
ray band (or both). We plan to re-analyse the data of these bursts, 
looking for evidence of either a break in the X-ray spectrum, or a 
very hard optical spectrum (this will depend on the assumed optical 
extinction). This will impact also on the assumed value of the Nn 
derived by fitting the X-ray data with a simple power-law model. 



In our scenario the observed optical and X-ray fluxes can be 
often ascribed to different processes. Therefore, when analysing the 
simultaneous optical to X-ray spectral energy distribution, some 
caution should be made in interpreting it as a single component 
connecting the optical and X-ray data. This has to be consistent 
with the analysis of the entire light curve in both bands, to be con- 
fident that they are produced by the same process. Only in such 
cases the host galaxy dust extinction can be compared to the X-ray 
absorption. Spectral information from NIR to UV would be crucial 
to this goal. 
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APPENDIX A: PHOMETRIC DATA REFERENCES 
References of the photometric data plotted in Figs.[T]|2j[3]|4j|5][6] 

- GRB 050318: Still et al. (2005); 

- GRB 050319: Wozniak et al. (2005), Mason et al. (2006), 
Quimby et al. (2006), Huang et al. (2007), Kamble et al. (2007); 

- GRB 050401: Rykoff et al. (2005), De Pasquale et al. (2005), 
Watson et al. (2006), Kamble et al. (2008); 

- GRB 050408: Foley et al. (2006), de Ugarte Postigo et al. 
(2007); 

- GRB 050416a: Holland et al. (2007), Soderberg et al. (2007); 

- GRB 050525a: Torii & BenDaniel (2005); Malesani et al. 
(2005), Chiang et al. (2005), Mirabal, Bonfield & Schawinski 
(2005), Homewood et al. (2005), Haislip et al. (2005), Green et 
al. (2005), Klotz et al. (2005a), Blustin et al. (2006); 

- GRB 050730: Sota et al. (2005), Holman, Garnavich & Stanek 

(2005) , Burenin et al. (2005), Klotz et al (2005b), D'Elia et al. 
(2005a), Bhatt & Sahu (2005), Kannappan et al. (2005), Pandey 
et al. (2006); 

- GRB 050801: Monard (2005), Fynbo et al. (2005c), Rykoff et 
al. (2006); 

- GRB 050802: Pavlenko et al. (2005), Fynbo et al. (2005e), 
Oates et al. (2007); 

- GRB 050820a: Cenko et al. (2006a); 

- GRB 050824: Sollerman et al. (2007); 

- GRB 050922c: Norris et al. (2005), Jakobsson et al. (2005b), 
Andreev & Pozanenko (2005), Durig & Price (2005), Henych et 
al. (2005), Novak (2005), Piranomonte et al. (2005), D'Elia et al. 
(2005b), Covino et al. (2005), Li et al. (2005); 

- GRB 051111: Butler et al. (2006), Yost et al. (2007); 

- GRB 060124: Romano et al. (2006), Misra et al. (2007); 

- GRB 060206: Wozniak et al. (2006), Stanek et al. (2007), Cur- 
ran et al. (2007a); 

- GRB 060210: Stanek et al. (2007), Curran et al. (2007b), 
Cenko et al. (2008); 

- GRB 060418: Melandri et al. (2006a), Cobb (2006a), Jeh'nek 
Kubanek & Prouza(2006), Koppelman (2006), Chen et al. (2006), 
Hafizov et al. (2006), Karimov (2006), Molinari et al. (2007); 

- GRB 060512: Mundell & Steele (2006), Cenko (2006a), 
Milne (2006), De Pasquale & Cummings (2006), Cenko & Baum- 
gartner (2006), Sharapov Djupvik & Pozanenko (2006); 

- GRB 060526: Campana et al. (2006b), French & Jelinek 

(2006) , Covino et al. (2006), Lin et al. (2006), Brown et al. (2006), 
Khamitov et al. (2006a), Morgan & Dai (2006), Khamitov et 
al. (2006b), Rumyantsev & Pozanenko (2006), Kann & Hoegner 
(2006), Khamitov et al. (2006c), Baliyan et al. (2006), Khamitov et 
al. (2006d), Khamitov et al. (2006e), Terra et al. (2006), Khamitov 
et al. (2006f), Rumyantsev et al. (2006), Dai et al. (2007), Thone et 
al. (2008); 

- GRB 060614: French et al. (2006), Schmidt Peterson & Lewis 
(2006), Cobb et al. (2006), Fynbo et al. (2006b), Delia Valle et al., 
(2006), Gal- Yam et al. (2006), Mangano et al. (2007b); 

- GRB 060729: Grupe et al. (2007); 



- GRB 060904b: Skvarc (2006), Oates & Grupe (2006), Mesch- 
eryakov et al. (2006), Cobb & Bailyn (2006), Greco et al. (2006), 
Soyano Mito & Urata (2006), Huang et al. (2006), Asfandyarov 
Ibrahimov & Pozanenko (2006), Klotz et al. (2008); 

- GRB 060908: Nysewander et al. (2006), Antonelli et al. 
(2006), Morgan et al. (2006), Cenko et al. (2008); 

- GRB 060927: Guidorzi et al. (2006), Torii (2006a), Ruiz- 
Velasco et al. (2007); 

- GRB 061007: Mundell et al. (2007); 

- GRB 061121: Page et al. (2006), Melandri et al. (2006b), 
Uemura Arai & Uehara (2006), Marshall Holland & Page 
(2006), Halpern Mirabal & Armstrong (2006a), Cenko (2006b), 
Torii (2006b), Halpern Mirabal & Armstrong (2006b), Efimov 
Rumyantsev & Pozanenko (2006a), Halpern& Armstrong (2006a), 
Halpern& Armstrong (2006b), Efimov Rumyantsev & Pozanenko 
(2006b), Cobb (2006b); 

- GRB 061 126: Perley et al. (2008a), Gomboc et al. (2008); 

- GRB 0701 10: Malesani et al. (2007), Troja et al. (2007); 

- GRB 070125: Cenko & Fox (2007), Xing et al. (2007), Ue- 
mura Arai & Uehara (2007), Greco et al. (2007), Yoshida Yanag- 
isawa & Kawai (2007), Terra et al. (2007), Mirabal Halpern & 
Thorstensen (2007), Updike et al. (2008), Chandra et al. (2008); 

- GRB 071003: Perley et al. (2008b), Cenko et al. (2008); 

- GRB 071003: Covino et al. (2008a); Cenko et al. (2008); 

- GRB 080310: Milne & Williams (2008a), Covino et al. 
(2008b), Chen et al. (2008), Garnavich Prieto & Pogge (2008a), 
Yoshida et al. (2008), Kinugasa (2008), Garnavich Prieto & Pogge 
(2008b), Urata et al. (2008a), Wegner et al. (2008), Hill et al. 
(2008), Cenko et al. (2008); 

- GRB 080319b: Li et al. (2008a), Milne & Williams (2008b), 
Urata et al. (2008b), Li et al. (2008b), Cwiok et al. (2008), Covino 
et al. (2008c), Wozniak et al. (2008), Swan Yuan & Rujopakarn 
(2008), Jelinek et al. (2008), Novak (2008), Krugly Slyusarev & 
Pozanenko (2008), Perley & Bloom (2008b), Tanvir et al. (2008), 
Bloom et al. (2008). 
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